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A Glasma flux-tube model has been proposed to explain strong elongation on pseudorapidity rj 
of the same-side 2D peak in minimum-bias angular correlations from ^Jsnn = 200 GeV Au-Au 
collisions. The same-side peak or "soft ridge" is said to arise from coupling of flux tubes to radial 
flow. Gluons radiated transversely from flux tubes are boosted by radial flow to form a narrow 
structure or ridge on azimuth. In this study we test the conjecture by comparing predictions for 
particle production, spectra and correlations from the Glasma model and conventional fragmentation 
processes with measurements. We conclude that the Glasma model is contradicted by measured 
hadron yields, spectra and correlations, whereas a two-component model of hadron production, 
including minimum-bias parton fragmentation, provides a quantitative description of most data, 
although ?7 elongation remains unexplained. 

PACS numbers: 12.38.Qk, 13.87.Fh, 25.75.Ag, 25.75.Bh, 25.75.Ld, 25.75.Nq 
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I. INTRODUCTION 

The systematics of measured hadron production and 
multihadron correlations in p-p and more-peripheral 
A-A collisions are described quantitatively by a two- 
component model including soft and hard components [1, 
By hypothesis the soft component arises from longi- 
tudinal fragmentation of participant projectile nucleons 
as a result of soft momentum transfers (e.g. diffractive 
scattering). The hard component arises from minimum- 
bias large-angle parton scattering and transverse frag- 
mentation as a result of (semi) hard momentum transfers. 
Those interpretations have been tested and elaborated in 
several studies [3-9]. 

Conventional theoretical descriptions of soft and 
hard components combine parton distribution functions 
(PDFs) [hard component, perturbative QCD (pQCD)] 
with "limitingfragmentation" (parton splitting cascade, 
DGLAP [To', |Tl|) at larger momentum fraction x and a 
phenomenological nonperturbative approach (e.g. Lund 
string model) at smaller x [l2|. Examples of Monte Carlo 
models combining both aspects are PYTHIA for p-p col- 
lisions [13 and HIJING for A-A collisions [14]. 

An alternative description of longitudinal particle pro- 
duction in more-central A-A collisions is based on the 
Color Glass Condensate (CGC) [15]. The CGC model 
invokes a statistical ensemble of classical color charges 
(sources) at larger x and a radiated classical color field 
at smaller x described as a Glasma. There is obviously 
a correspondence between Glasma vs Lund strings at 
smaller x, and color radiators vs pQCD parton splitting 
cascade at larger x. We wish to explore those dichotomies 
by quantitative comparisons of theory with spectrum and 
two-hadron correlation data. 

In a previous study we considered the relation between 
Glasma flux-tube predictions and angular-correlation 



phenomenology for same-side (defined below) correla- 
tions in p-p and Au-Au collisions, specifically the so- 
called "soft ridge" 0. The study concluded that 
whereas pQCD-based descriptions of jet-related angular 
correlations are in quantitative agreement with data the 
description based on Glasma flux tubes is inconsistent 
with correlation data in several ways. 

In the present study we consider the relation between 
Glasma predictions [l3 measured two-particle 
(gluon and hadron) correlations on transverse momen- 
tum pt (or transverse rapidity yt defined below) as well 
as pseudorapidity r] and azimuth (j). We find further sub- 
stantial discrepancies between CGC-Glasma predictions 
and data. We conclude that even if conjectured radial 
flow played a role in nuclear collisions the Glasma flux- 
tube model does not correspond to physical processes 
that might, in the presence of such flow, produce the 
observed same-side 2D peak structure elongated on rj. 

This article is arranged as follows: Sec. HIl re views anal- 
ysis methods applied to RHIC hadron data. Sec. [TTIl de- 
scribes gluon correlations predicted by the Glasma flux- 
tube model. Sec. IIVI describes measured hadron correla- 
tion and spectrum data, and Sec. |V] compares predicted 
Glasma gluon correlations with the hadron data. 



II. ANALYSIS METHOD 

We review technical aspects of STAR correlation anal- 
ysis applied to nuclear collisions at the Relativistic Heavy 
Ion Collider (RHIC). Method details are provided in 
Refs. [i&Si, [19^23. Topics include A-A collision ge- 
ometry, the two-component model of hadron production, 
correlation measures, 2D histograms, model functions for 
2D fits and the relation between fluctuations and an- 
gular correlations. 
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A. A-A collision geometry 

A-A collision geometry is described by the Glauber 
model relating the A-A differential cross section to par- 
ticipant nucleon number N^art ctnd N-N binary-collision 
number A^^^^ [23|. A derived projectile- nucleon mean 
path length v = 2Ni)in/Npart is also defined. Through 
the measured A-A differential cross section on charged- 
hadron multiplicity rich within some angular acceptance 
the Glauber parameters can be related to observed rich- 

Optical Copt [21] and Monte Carlo cmc [24] eccen- 
tricities have been invoked to model A-A eccentricity 
required for interpretation of the azimuth quadrupole 
measured by V2. The former assumes a smooth mat- 
ter distribution across nuclei whereas the latter assumes 
that point-like participant nucleons are the determining 
elements. A priori support for Copt assumes that the 
azimuth quadrupole emerges from interactions at small 
X < 0.01 where one might expect onset of a smooth, 
saturated glue system (e.g. Glasma) [25]. A posteriori 
support for Copt arises from a simple systematic trend 
(X Nijinclp^ observed for V2{2D} data which accurately 
exclude contributions from jet structure (nonflow) [21]. 

B. Two-component hadron production model 

According to the two-component model spectra and 
correlations from nuclear collisions can be decomposed 
(near mid-rapidity) into soft and hard components de- 
noting respectively longitudinal fragmentation (mainly 
diffr active dissociation) of projectile nucleons and trans- 
verse fragmentation of large- angle-scattered partons [l|,[3|, 
0, 0] • Soft and hard components from p-p collisions are 
clearly distinguishable in pt x pt or yt x yt correlations 
(defined below). 

In more-peripheral A-A collisions the soft component 
should vary with centrality ex Npart and the hard compo- 
nent should vary ex Ni^in [26], those trends constituting 
the Glauber linear superposition reference for spectra and 
correlations. The soft-component correlation amplitude 
drops to zero in high-multiplicity p-p and more-central 
Au-Au collisions. The hard-component fraction is a few 
percent in minimum-bias (NSD) p-p collisions [f] but in- 
creases to about one third of the final-state hadron yield 
in central Au-Au collisions [5]. Thus, the two-component 
model of hadron spectra and correlations provides a self- 
consistent quantitative description of almost all aspects 
of nuclear collisions based on pQCD and measured prop- 
erties of elementary e+-e~ and p-p collisions. 



C. Correlation measures 

Two-particle correlations are structures in pair- 
density distributions on six-dimensional momentum 
space {pti , ^1 , 01 , Pt2 , ^2 , 02 ) • We visualize correlation 
structure in 2D subspaces (pt^Pt) and (t^a^^a) (defined 



below) which retain almost all structure within a lim- 
ited T] acceptance such as the STAR Time Projection 
Chamber (TPC) [13 • We measure correlations with 
per-partide statistic Ap/y^p^ = Po ((^) — 1), where 
Ap = p — p^^f is the correlated-pair density, pref is 
the reference- or mixed-pair density, (f) is the (unit- 
normal) sibling/mixed pair number ratio and prefactor 
Po = nch/Ar]A(l) is the charged-particle 2D angular den- 
sity averaged over angular acceptance {Arj^Acj)) [ol, |28|. 
Pair ratio f is averaged over kinematic bins (e.g. multi- 
plicity, Pt, vertex position), and we assume factorization 
of the reference pair density pref ^ Po- 

The per-particle measure is an improvement over con- 
ventional per-pair correlation function (f) p/ pref or 
(f) — 1 Ap/ pref since it eliminates a trivial l/uch trend 
common to all per-pair measures (except for quantum 
correlations). The intensive definition in terms of hadron 
2D density po(^) (^ is the A-A impact parameter) rather 
than multiplicity rich also eliminates a trivial dependence 
on detector angular-acceptance factor ArjAcj). |2Qj . 



D. Transverse- momentum correlations on pt x pt 

2D correlations on pt or transverse rapidity yt = 
ln[(pt + mt)/m^] {m^ is assumed for unidentified 
hadrons) are complementary to 4D angular correlations 
in 6D two-particle momentum space, yt is preferred for 
visualizing correlation structure on transverse momen- 
tum. Pt X Pt or yt X yt and angular correlations can 
be defined for like-sign (LS) and unlike-sign (US) charge 
combinations and also for same-side (SS) and away-side 
(AS) azimuth subregions of angular correlations (defined 
below). Manifestations of different correlation mecha- 
nisms (e.g. soft and hard components) can be clearly dis- 
tinguished in the four combinations of charge-pair type 
and azimuth subspace, with distinctive forms for each of 
the LS and US charge combinations and for SS and AS 
azimuth subspaces 0, 0] • 

For correlations on pt x pt the pair ratio is r{pti^Pt2) 
and the prefactor becomes ^/p^ = \/ Po{Pti) Po{Pt2)^ the 
geometric mean of single-particle pt spectra. The correct 
prefactor is essential for proper comparison of predicted 
gluon correlations and measured hadron correlations. 

E. Angular correlations on (?7a,0a) 

Angular correlations can be formed by integrating over 
the entire pt x pt pair acceptance (minimum-bias angular 
correlations) or over subregions 0, • Examples of the 
latter include "trigger-associated" dihadron correlations 
resulting from asymmetric cuts on pt x pt ^] . 

Two-particle angular correlations are defined on 4D 
momentum subspace (t^i, 7^2, ^i, ^2)- Within acceptance 
intervals where correlation structure is invariant on mean 
angle (e.g. rjj^ = ^1 +^2) angular correlations can be pro- 
jected by averaging onto difference variables (e.g. t^a = 
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Vi ~ V2) without loss of information to form angular au- 
tocorrelations [20\^2§\. The 2D subspace (77a7(/>a) is then 
visuahzed. Symbol Ax is used as a measure of the detec- 
tor acceptance on parameter x. 

Angular correlations can be formed separately for like- 
sign and unlike-sign charge combinations, as well as 
for the charge-independent (CI = LS + US) combina- 
tion [9, 28]. The pair angular acceptance on azimuth can 
be separated into a same-side (SS) region (|^a| < 7r/2) 
and an away-side (AS) region (|0a| > 7r/2). The SS re- 
gion includes intranet correlations (hadron pairs within 
single jets), while the AS region includes inter }ei corre- 
lations (hadron pairs from back-to-back jet pairs). 



particle sources. In the NBD context Acr^(fe, Ax) = n/k 
(correlated particles per independent source) may then 
represent the integral up to a specific bin size (scale) of 
angular correlations arising from a superposition of sev- 
eral physical mechanisms within a specific acceptance. 
n/k typically increases monotonically with increasing bin 
size, angular acceptance and correlation amplitude. Sys- 
tematic details may refiect several underlying correlation 
mechanisms. Only differential correlation analysis can 
identify individual correlation sources (see Sec. lIVp . 

III. GLUON CORRELATIONS FROM GLASMA 



F. Angular-correlation model functions 

The hard component of angular correlations includes a 
SS 2D peak at the angular origin and an AS ID peak on 
azimuth uniform on ?7a (within the STAR TPC accep- 
tance). The minimum-bias SS 2D peak (intrajet correla- 
tions) is well modeled by a 2D Gaussian. Except for p-p 
and more-peripheral A-A collisions the AS peak (inter- 
jet correlations) is conveniently modeled as an AS dipole 
(X cos(^A — tt). The soft component is modeled by a ID 
Gaussian on ?7a with r.m.s. width ^ 1. 

The combined model function including azimuth 
quadrupole term cos(2^a) required to describe A-A an- 
gular correlations within the STAR TPC is [d, §, ^ 

= Ao+AiDe"H^)' (1) 

y/Pref 

+ A2De"H(^) ^^^^ } 

+ Ad [1 + cos(0A ~ tt)] /2 + Aq 2 cos(2 ^a)? 

where a narrow 2D exponential describing quantum cor- 
relations and electron pairs from 7 conversions [6] has 
been omitted for clarity. The nonjet quadrupole ampli- 
tude is expressed in terms of conventional parameter V2 
by Aq{2D} = po{b)vU2D}{b) M- 

G. Fluctuations from correlations 

The direct relation between fluctuations and two- 
particle angular correlations was established in Ref. (2Q| 
and implemented for (pt) fluctuations in [30|, [3l|. Given 
variance = (n^) — the per-particle number fluctu- 
ation measure A(j^(fe, Ax) = (a^ — n)/n is equivalent 
to the integral up to scale (bin size) 5x within some ac- 
ceptance Ax of per-particle number angular correlation 
measure Ap/ ^J~p^ deflned on difference variable xa [20]. 

The negative binomial distribution (NBD) is a two- 
parameter multiplicity distribution for some bin size bx 
over some acceptance Ax (which should be specifled). 
NBD parameters are mean multiplicity fi and parameter 
k which can be interpreted as a number of independent 



Description of so-called "bulk" hadron production in 
terms of a CGC Glasma is an alternative to conven- 
tional parton distribution functions (PDFs) and longi- 
tudinal fragmentation of projectile nucleons. According 
to Ref. [17] bulk hadron production arises from low-x glu- 
ons residing in the wavefunctions of projectile nucleons, 
and "leading correlations are present already in the wave 
functions of the colliding objects." 

If a characteristic energy scale (saturation scale Qs) 
is large enough (Qs ^ ^qcd) the longitudinal system 
may be divided into radiating color charges at larger 
X (sources) and a (nominally boost invariant) saturated 
classical fleld (Glasma) at smaller x. The transverse cor- 
relation length for the Glasma fleld is 1/Qs, interpreted 
as the transverse size of Glasma "flux tubes." 

Some features of angular correlations in more-central 
Au-Au collisions, particularly ri elongation of the SS 2D 
peak, have been interpreted as a "signal" for Glasma flux 
tubes in the A-A initial state (IS) [ITI, |32|. The present 
analysis is intended to test the conjecture that ri elonga- 
tion of the same-side 2D peak in more-central Au-Au col- 
lisions can be explained in terms of interaction of initial- 
state Glasma flux tubes with radial flow. 



A. Glasma gluon fluctuations and the NBD 

In the Glasma model NpT = Q\^^ ^ 1 is the 
number of flux tubes in an A-A system with trans- 
verse area [17], and NpT = 0(1) A/p^rt [IS]. Each 
flux tube emits into — 1 color states for a total of 
(N^ — 1)Nft/2ti independent Bose-Einstein color radia- 
tors. Each radiator emits on average fig = 0(1) / as{Ql) 
gluons. In Eq. (3.2) of Ref. [17] the gluon total multiplic- 
ity is given by Ng = {/n /as)QlS±. 

The Glasma model is thus a two-tiered statistical sys- 
tem. Fluctuations in the Glasma gluon multiplicity can 
be modeled by a negative binomial distribution with pa- 
rameters Ng and k. According to the Glasma model in 
Ref. [17] multigluon correlations in the classical limit ex- 
pressed in terms of the variance of gluon number Ng 
have Aa%^ = {a%^ — Ng)/Ng = Ng/k = or no cor- 
relations. Averaging over A-A color source conflgura- 
tions induces signiflcant correlations, with Ng/k ^ 0. 
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Correlations in p-p collisions are expected to be small 
{k oo). However, that trend is inconsistent with ob- 
servations (Sec. IIV A[) . Glasma gluon fluctuations are 
further discussed in Sec. IVI Gi 



B. Glasma predicted two-gluon correlations 

Perturbative results from Ref. [18] are summarized by 

C2{Pl^P2) Ap{pi,p2) = p{pi,P2) - Pref{Pl^P2) 

Ap{Pl^P2) _ 1 _ 1^2 
Pref{pi,P2) k QlS^ 

Pref{pi,P2) P0{P1)P0{P2) 

poip) (X Hpt/Qs){Qs/Pt)\ (2) 

with NBD parameter k ^ {N^ - 1)Nft/27t flT^ and 
1^2 = NpT/k ^ 27r/(7V^^ - 1) = 0(1) approximately 
constant for pt/Qs ^ 1 [18]. Since k ~ Nft we have 
A^^ ^ figNpT- The Glasma predicted energy dependence 
is /c ~ Q? ~ Js [17]. As defined, hZ2 has the per- 
particle [gl [ol, [20| structure nc/i(^ — 1) ^ Nft{^ ~ 
For the perturbative case the only correlation source is 
production hierarchy: multiple gluons associated with 
each independent fiux-tube radiator. Perturbative two- 
gluon correlations are then independent of angle differ- 
ences and factorizable on {pti^Pt2)- A non-perturbative 
calculation is required to determine correlations down to 
small pt [18]. hZ2 may then depend on pt and relative 
angle 0a , but not on tja if fiux tubes are boost invariant. 

IV. MEASURED HADRON CORRELATIONS 

There is now an extensive phenomenology of hadron 
correlations from which the underlying parton dynam- 
ics may be inferred. The phenomenology described 
here exhausts all differential correlation structure in p-p 
and Au-Au collisions, and therefore all fiuctuation phe- 
nomenology. Fluctuation measures, as running integrals 
of differential correlations, are consistent with correlation 
measurements but retain less information. 



A. Correlations from p-p collisions 

Figure [1] shows parametrizations of yt x yt and 2D an- 
gular correlations from p-p collisions at 200 GeV [1, Q . 
Charged-hadron correlations on yt x yt (left panel) with 
yt G [1,4.5] or pt G [0.15,6] GeV/c integrated over the 
STAR TPC angular acceptance \r]\ < 1 and 27r azimuth 
include two well-separated peaked structures described 
as the soft component and the hard component Q. 
The soft component is a 2D peak localized below 0.5 
GeV/c (e.g. within yti + yt2 < 2). The hard component 
is a second distinct 2D peak which dominates the comple- 
mentary yt X yt space. The peak is centered near yt ^ 2.7 



{pt ~ 1 GeV/c) and does not extend below 0.35 GeV/c in 
p-p collisions. Corresponding structures are superposed 
in minimum-bias angular correlations (right panel). 




FIG. 1: (Color online) Two-particle hadron correlation 
histograms from 200 GeV p-p collisions based on data 
parametrizations from Refs. 0|. Left: yt x yt correlations 
within the angular acceptance showing soft and hard com- 
ponents, Right: pt-integral angular correlations on (?7a,(/>a) 
showing a superposition of soft and hard components. 

Figure [2] shows distinct structures in angular corre- 
lations on (7^a7</>a) corresponding to the soft and hard 
components separated by cuts on pt x pt or yt x yt. 
The soft component of angular correlations (left panel) 
is a ID Gaussian on tja with approximately unit r.m.s. 
width characteristic of "short-range" correlations. The 
soft component is almost entirely US pairs. Such struc- 
ture is consistent with diffractive scattering and longi- 
tudinal fragmentation of projectile nucleons to hadrons 
(near mid-rapidity) in charge-neutral (US) combinations. 




FIG. 2: (Color online) Two-particle hadron correlation 
histograms from 200 GeV p-p collisions based on data 
parametrizations from Refs. Q, Left: Soft component 
(longitudinal projectile nucleon fragmentation to US pairs), 
Right: Hard component (large- angle- scattered parton frag- 
mentation to US SS pairs and CI=LS+US AS pairs). 

The hard-component peak on {yt^yt) corresponds to 
two structures on (t^a^^a) (right panel), a SS 2D peak 
centered at the angular difference origin and an AS ridge 
centered at tt on ^a and uniform on 77a- The SS 2D peak 
is almost entirely US pairs (reflecting local charge conser- 
vation during fragmentation), whereas the AS ID peak 
on azimuth (jet-jet ridge) is composed of equal numbers 
of LS and US pairs (no interjet charge correlation). 

The structure on (yt^yt) corresponding to the SS 2D 
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peak lies close to the main diagonal {yn ^ yt2)^ whereas 
that corresponding to the AS ridge is substantially broad- 
ened relative to the main diagonal [1, • Those system- 
atics are consistent with large-angle parton (gluon) scat- 
tering and fragmentation to back-to-back jet pairs, with 
hadron local charge and momentum conservation. 

B. Correlation evolution with A- A centrality 

If A-A collisions were simply linear superpositions of 
p-p (or N-N) colhsions (eikonal model) we could extrap- 
olate p-p correlation phenomenology according to the 
Glauber model to describe A-A data. The eikonal model 
provides a reference system denoted by Glauber linear 
superposition (GLS), including participant scaling of the 
soft component and binary-collision scaling of the hard 
component [s*, "26^. Correlation data from Au-Au colli- 
sions at 62 and 200 GeV follow the GLS reference from 
peripheral collisions to an intermediate centrality and 
then transition to substantially different behavior 
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FIG. 3: Left: Amplitude of the same-side 2D Gaussian fitted 
to minimum-bias 2D angular correlation data from 200 GeV 
Au-Au collisions 6]. Right: Fitted peak widths for the same- 
side 2D Gaussian. GLS indicates a Glauber linear superposi- 
tion reference extrapolated from measured p-p collisions 1]. 

Figure [3] summarizes centrality variation of SS 2D peak 
parameters. Deviations from the GLS reference extrapo- 
lation above the transition point include a rapid change 
(within one 10% centrality bin) in the rate of increase 
(slope on centrality measure u defined in Sec. Ill A[) of SS 
and AS jet-related peak amplitudes, a rapid increase in 
the SS 2D peak t^a width and a significant decrease in 
the width. The SS 2D peak aspect ratio transitions 
from nearly 2:1 elongation on azimuth to 3:1 elongation 
on pseudorapidity Q . Jets in p-p collisions are nearly as 
anomalous as those in central Au-Au collisions. Those 
trends have been interpreted quantitatively in terms of 
modified parton fragmentation in more-central Au-Au 
collisions [8]. 

Figure H] shows jet structure in 200 GeV central {b = 0) 
Au-Au collisions extrapolated from measured centrality 
trends. The nonjet quadrupole measured by V2i2D} 
is zero with small upper limit in that case [21|. The 
AS ID peak on azimuth is uniform on tja within the 
STAR TPC acceptance, as is the nonjet quadrupole. In 




FIG. 4: (Color online) Left: Angular correlations from 200 
GeV Au-Au collisions with centrality ly — 6 {b — 0) with 
nonjet fit components (quadrupole, ID Gaussian on ?7a) sub- 
tracted to reveal nominal jet correlations. The vertical zero 
is the estimated true zero offset for these overlapping jet cor- 
relations. Right: The previous histogram with the away-side 
dipole term subtracted to isolate the same-side 2D jet peak 
extrapolated to 77 G [—2, 2]. 

Fig. |4] (right panel) the r^A-independent components have 
been subtracted, leaving the SS 2D peak (extrapolated to 
|^a| =4) as the remaining correlation component. The 
SS 2D peak for minimum- bias (p^-integral) data is always 
consistent with a 2D Gaussian (no additional structure). 
Although "trigger-associated" dihadron correlations do 
reveal non-Gaussian (on t^a) SS features, those data are 
subsets of the minimum-bias ensemble and should reflect 
the same basic fragmentation process. 



C. Relation to hadron spectra and yields 

Soft and hard components of two-particle correlations 
have counterparts in single-particle hadron pt spectra. 
The p-p Pt spectrum hard component [l[ corresponds 
quantitatively to the hard component in Fig. [1] (left 
panel), p-p spectra and correlations form a simple system 
quantitatively consistent with all aspects of pQCD down 
to zero hadron momentum and 3 GeV parton energy [H . 

In a minijet context the pt spectrum hard component 
(single-particle fragment distribution) is the marginal 
projection of the correlation hard coniponent on pt x pt 
or yt X yt (fragment pair distribution [^). For all Au-Au 
centralities the measured correlation hard component on 
PtXPt [33] is quantitatively consistent with the hard com- 
ponent inferred from measured pt or yt spectra [3I] . The 
hard-component peak on pt xpt or yt x yt persists as a dis- 
tinct structure with mode near pt = I GeV/c even in cen- 
tral Au-Au collisions. For both spectra and correlations, 
suppression at larger pt ^ 10 GeV/c is accompanied by 
much larger enhancement at smaller pt ~ 0.5 GeV/c Q. 
Suppression and enhancement trends on Au-Au central- 
ity are closely (anti) correlated j8|. 

The integral of the SS 2D peak in angular corre- 
lations (all jet-related hadron pairs) combined with a 
pQCD dijet cross section can be converted to a hard- 
component hadron yield. Variation of the calculated 
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hard- component yield with centrahty explains evolution 
of the Au-Au total hadron yield with centrality, revealing 
that about one third of the total yield in central Au-Au 
collisions is included in resolved minijets [5]. Analysis of 
the pt or yt dependence of the SS 2D 7^-elongated peak 
in angular correlations [H, [H, [33| shows that this struc- 
ture corresponds to the hard-component peak on pt x pt 
or yt X yt near pt = I GeV/c for all Au-Au centralities. 



V. GLASMA FLUX TUBES vs HADRON DATA 

The Glasma model is a one-component (soft) model of 
gluon production near mid-rapidity emphasizing more- 
central A-A collisions. There is no (semi)hard parton 
(gluon) scattering to mid-rapidity. The Glasma model 
competes with the two-component (soft plus hard) model 
of hadron production inferred from yields, spectra and 
two-particle correlations derived from p-p and A-A col- 
lisions E, S i, B S i- We provide direct comparisons 
between Glasma predictions and measured hadron spec- 
trum and correlation data as a test of Glasma relevance 
to the hadronic final state in nuclear collisions. We con- 
sider spectrum and correlation trends on (7^, ^), pt or yt^ 
A-A centrality and collision energy. 




(GeV/c) 



2 4 6 8 10 
Pj,mj-m^ (GeV/c) 



FIG. 5: Left: Gluon pt spectrum from Glasma flux- tube 
model (with added factor p?). Right: Comparison of hadron 
mt spectrum soft component from Ref. ^ (dashed curve) and 
gluon Pt spectrum from Glasma flux-tube model (solid curve) . 



T] density for the shape comparison. The hard gluon spec- 
trum has no observed A/p^r^-scaling counterpart in the 
hadronic final state. There is no correspondence in the 
gluon spectrum to the measured hard component scaling 
as Ni)in in hadron spectra [sj or to the large-angle- 
scattered parton spectrum varying as 1/Pt'^ [Si]- 



B. Glasma and NBD k-parameter systematics 



A. Glasma single-gluon spectrum 

The number of Glasma flux tubes is expected to scale 
with A-A centrality as QlS± ~ Npart [16]. The cor- 
responding gluon production should then be compared 
with the soft component of hadron production which 
does follow Npart scaling ^ Q . The inclusive single- 
gluon spectrum from Ref. [18|, Fig. 3 (left panel) is well 
represented by 

^ = erf(p,/0.3 GeV/c) r, (3) 

dp't ' ^ {(1/5)2 + (p?/33)2}^'^ ^ 

whereas the soft-component charged-hadron spectrum 
for p-p and all centralities of A-A is the Levy distribution 

2 d^Uch ^ 45 

ripart drjmtdmt {1 + (mt - mT,)/nT}'^' 

with T = 0.145 GeV and n = 12.8 [l|,ll. 

Figure [5] (left panel) shows Eq. (j3j) which compares 
well with Fig. 3 (left panel) from Ref. [l8| where the 
spectrum has been scaled by factor p^ to show details 
at smaller pt. The Glasma expectation for the gluon 
spectrum trend is ex \n{pt/Qs){Qs/pt)^ at larger pt/Qs 
with Qs ^ 1 GeV [18]. 

Figure [5] (right panel) shows a comparison between the 
Glasma gluon spectrum ex 1/Pt and the per-participant- 
pair hadron spectrum soft component cx l/m}'^-^ at larger 
Pt and rrit respectively. The mt spectrum integrates to 
duch/dr] = 2.5 (NSD p-p collisions). The gluon spectrum 
in Eq. (j3|) has been divided by 6 to approximate the same 



In Fig. [6] (left panel) we show k data from Fig. 3 
of Ref. [17] plotted as l//c, which then increase by fac- 
tor 2 from RHIC to LHC energies. For "short-range" 
(localized on r]) angular correlations we do expect 1/k 
(representing per-pair angular correlations) to decrease 
with increasing system size (e.g. A-A centrality). In 
the CGC model the k ^ parameter is expected to 
increase monotonically with collision energy, or equiv- 
alently 1/k should decrease with collision energy. The 
observed strong 1/k increase with energy contradicts the 
Glasma expectation. 




10 10^ 10^ 10 10 10^ 10^ 10 



Vs (GeV) Vs (GeV) 

FIG. 6: Left: Energy dependence of negative binomial dis- 
tribution (NBD) parameter k from Ref. 17] plotted as 1/k 
for p-p collisions and angular acceptance |?7| < 1 or Ar] = 2, 
Right: The same data plotted as drich/dr] xl/k measuring the 
integral of angular correlations within the anglar acceptance. 
The curves are based on In (y^/ 13. 5 GeV) (see text). 

In Fig. [6] (right panel) corresponding product 
drich /df] X 1/k measures per-particle angular correlations 
integrated on scale up to the angular acceptance (in this 



7 



case Af] = 2 and = 27r). We conclude that integrated 
angular correlations in p-p collisions from whatever mech- 
anism increase with collision energy faster than \og{^/s). 

Given the measured energy dependence of pt angu- 
lar correlations [30], nonjet azimuth quadrupole cor- 
relations [2l| and minijet angular correlations [6] in 
p-p and Au-Au collisions at and below 200 GeV we 
characterize the k data based on the energy trend 
In (y^/ 13. 5 GeV). The solid curve in the right panel 
is drich/dr] x 1/k = 0.18{ln(yi/13.5 GeV)}^-^^ The 
hadron yield increase with energy is well described 
above 200 GeV by duch/dr] ^ 0.88{ln(v^/13.5 GeV)}. 
The solid curve in the left panel is therefore 1/k = 
0.21{ln(v^/13.5 GeV)}°•'^^ generally consistent with 
QCD processes (e.g. minijet production) but inconsistent 
with the flux-tube expectation k oc ^/s^ with A > [17\ . 

We can express the energy and centrality dependence 
of Glasma correlation parameter in terms of mea- 
sured quantities. Since i<i2 = NpT/k ~ Npart/k we 
have ^ {Npart/'^nch) X {"^nch/k). Combining the 
two-component expression for the flrst factor [26] with 
the trend inferred from Fig. [6] (right panel) we obtain 
^ 0.36{ln(v^/13.5 GeV)}i '^V{l + O-U^ - !)}• sub- 
stantial increase with energy and decrease with centrality. 



C. Glasma gluon pt x pt correlations 

Figure [7] (left panel) shows Fig. 2 (right panel, AS 
pairs) of Ref. [17]. Fig. 2 (left panel, SS pairs) of Ref. [I^ 
includes self pairs along the diagonal but is otherwise sta- 
tistically equivalent to the right panel. The cited flgure 
shows a plan view of i<i2 on {pti-,Pt2) relative to saturation 
scale Qs ^ 1 GeV/c. The histogram is replotted here in 
isometric view on pt in units GeV/c. 

Fig. [7] (left panel) is analogous to correlations plotted 
on yt X yt (Fig. [H left panel). K2 prefactor QlS±_ in Eq. 
(5.1) of Ref. is equivalent to rich in ^ch{'^{Pti-,Pt2) — l} 
in the notation of Sec. Ill CI Thus, K.2 should be com- 
pared with per-pair measure Ap/ prefiyt^Vt) which 
increases dramatically at larger yt^ strongly contradict- 
ing the Glasma prediction plotted in Fig. [71 (left panel). 




FIG. 7: (Color online) Left: Data from Fig. 2 (right panel) of 
Ref. [13 plotted in isometric view showing predicted Glasma 
t wo-gluon corre lations. Right: The same data with prefactor 
^/pref{pti,Pt2) based on a hadron spectrum soft component. 



Figure [71 (right panel) shows the k,2 histogram in 
the left panel multiplied by prefactor \/pref{PtiiPt2) = 

\/ po{pti)po{pt2)^ where po{pt) is the hadron soft compo- 
nent in Eq. @ divided by 27r to form a 3D density. Pref- 
actor po{b) = rich/ ArjAcj) appropriate for angular correla- 
tions is replaced by the geometric mean of single-particle 
Pt or yt spectra appropriate for pt x pt correlations. The 
plotted histogram is proportional to per-particle measure 
^pI \/~Pref{pti^Pt2) and thus directly comparable with 
per-particle correlation data. The result is consistent 
with measured soft-component pt x pt correlations. 

Figure [51 (left panel) shows the Glasma data in Fig. [71 
(left panel) with the prefactor formulated using the gluon 
spectrum deflned by Eq. (|3|), which is much harder than 
the hadron spectrum. The result is still generally consis- 
tent with soft-component pt x pt or yt x yt correlations. 




FIG. 8: (Color online) Left: Glasma predictions in the form 
Ap/pref from Fig. [71 (left panel) multiplied by prefactor 
^/pTIf = \/po(pti, b)po{pt2, b) (gluon spectra) to obtain the 
per-particle form Ap/y/pref- Right: The yt x yt histogram in 
Fig. [1] (left panel) transformed to pt with the proper Jacobian. 

Figure [51 (right panel) shows the p-p parametrization 
in Fig. [H (left panel) representing data in Refs. |7j re- 
plotted on {pti^Pt2) with Pt G [0.2,5] GeV/c and with 
the proper Jacobian. The hard-component peak mode is 
near 1 GeV/c as expected, and the soft-component peak 
is just visible near the origin. The shape of the hard com- 
ponent on Pt xpt projected to ID is in quantitative agree- 
ment with the hard component inferred directly from p-p 
single-particle pt spectra in Ref. [1] (Fig. 10, left panel). 

This comparison reveals that the Glasma flux-tube 
model is strongly contradicted by measured hadron pt xpt 
correlations. The Glasma histogram corresponds qualita- 
tively to the soft component of measured hadron correla- 
tions. There is no corresponding hard-component peak, 
no large-angle parton scattering, in the Glasma model. 
The qualitative difference persists even when a prefactor 
derived from a hard gluon spectrum is introduced. 

One might argue that the intervening hadronization 
process could invalidate such a comparison. The max- 
imum hadron momentum per gluon would result from 
direct 1^1 correspondence as in local parton-hadron 
duality (LPHD) [35] (but observed local charge conser- 
vation is not then respected). Fragmentation (1 ^ 2 or 
more) should actually reduce the mean hadron momen- 
tum relative to the Glasma spectrum. There is thus no 
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possibility, within the Glasma model of "bulk" hadron 
production, to generate a counterpart to the observed 
Pt X pt or yt X yt hard component. 

D. Glasma pt x pt correlations and radial boosts 

A proposed mechanism for formation of the SS 2D peak 
in angular correlations from more-central Au-Au colli- 
sions is radial boost of Glasma flux tubes. We can then 
ask what would be the effect of conjectured radial flow on 
the predicted i^2{pti-,Pt2) in Fig. [71 (left panel)? Does ra- 
dial flow also produce the nominal hard-component struc- 
ture in Pt X Pt correlations, as in Fig. [8] (right panel)? 

Starting with the 2D histogram in Fig. [3 (left panel) 
the following procedure was applied with two pt spec- 
trum models po{pt) (prefactor Nft is ignored to simplify 
terminology, 1^2 Ap/pref)- Reference pref(jPti,Pt2) = 
po{pti)po{pt2) was formed and numerator Ap = pref^2 
was then obtained from the histogram in Fig. [71 (left 
panel). Azimuth angles were randomly sampled. Ap 
and pref were boosted from {pti^Pt2) to (p^i,p^2) t)y 
{Pt) = 0-6 (maximum mean value inferred from cen- 
tral Au-Au collisions [36^). Boosted Ap/ p^ef 1^2 was 
then recovered. The actual procedure was based on a 
Monte Carlo sampling (10^ samples) of ^^2, {pti^Pt2) and 
(01, ^2) to construct boosted histograms. The two po{pt) 
spectrum shapes were a Maxwell-Boltzmann (M-B) with 
Teff = 0.15 GeV and the gluon spectrum of Eq. ([3]). 




FIG. 9: (Color online) Correlations from Fig. [71 (left panel) 
boosted by mean radial speed = 0.6 (see text). Left: 
Single-particle spectrum is Maxwell-Boltzmann with Teff = 
0.15. Right: Single-particle spectrum is defined by Eq. ((3|). 

Results are shown in Fig.[9lfor the two cases. For those 
boost conditions the hZ2{pti^Pt2) distribution is changed 
modestly, mainly an increase at larger pt from 1 to 2 
or 2.5 which cannot possibly match the hard-component 
structure evident in Fig. [51 (right panel) when converted 
to the per-particle form y/ Pref{pti,Pt2) i^2{pti,Pt2)' 

E. Glasma gluon azimuth correlations 

Figure [TOl (left panel) reproduces Fig. 5 (right panel) 
of Ref. [17] which is related indirectly to angular corre- 
lations plotted on (t^a^^a) 0, B S Isl- The left panel 



in Ref. \y?\ includes extraneous self pairs ("large delta 
function peaks") along one axis which are removed in 
the right panel. The horizontal axes are magnitudes of 
sum and difference vectors pti + pt2 and pti — pt2 rel- 
ative to Qs = 1 GeV. One axis corresponds to parallel 
pairs, the other to antiparallel pairs. Correlation struc- 
ture seems to indicate a preference for momenta parallel 
and antiparallel (near the axes). The plot is interpreted 
to indicate that a "collimation effect" might be present 
in the initial state of p-p collisions, which might in turn 
explain the "ridge" observed in 7 TeV p-p collisions [s^ • 




FIG. 10: (Color online) Left: Histogram from Fig. 5 (right 
panel) of Ref. fl7] showing two-gluon correlations on vector- 
momentum sum and difference. Right: The Jacobian from 
to cos((/)a) plotted as a surface on (x+,x-) (see text). 

The coordinate axes in Fig. [TOl represented by X- 
are defined by 

xl = \pti±pt2\^ = 4p?[l±sin(2V;)cos((/)A)]/2, (5) 

where pti — V^Pt sin{ijj) and pt2 = V^Pt cos(?/^) define pt 
and i/j (polar coordinates in Fig. [7|). The axes extend to 
'^Pt,max = 6 GeV/c, twice the maximum single-particle 
momentum in the pair sample. The plotting variables can 
be inverted to sin(2?/^) cos(^a) = {x\ — x'^_)/{x\ + x?_) 
(normalized difference diagonal) and Apl = x\^x'^_ (sum 
diagonal). The sum diagonal measures the quadratic 
mean of two transverse- momentum magnitudes. The dif- 
ference diagonal measures angle ^a between pairs of mo- 
menta. For projection to space (x+,x_) an average over 
sin(2?/^) would be determined by the (slowly varying) dis- 
tribution in Fig. [71 (left panel). There is no sensitivity to 
elongation on t^a, to a SS ridge per se. 

The apparent correlation structure is symmetric about 
(/>A = 7r/2, and the angular distribution seems to he 
sharply peaked toward the AS limit as well as toward 
the SS limit ("collimation effect even in the absence of 
radial flow" [17]), but the plot is misleading. Figure 9 
(right panel) of Ref. [18] shows the same correlation mea- 
sure i<i2 plotted directly on (j)/^ for pti ~ Pt2 ^ 3 GeV/c 
[therefore sin(2?/;) ^ 1]. The distribution on ^a is nearly 
uniform (variation within ±10%). The relation between 
plots on (x+,x_) and on (j)/^ appears to be the Jacobian 
1/ sin(0A) from (j)/^ to cos(^a)- 

Figure [TOl (right panel) shows the surface l/sin(^A) 
plotted on (x+,x_). The detailed agreement with the 
left panel is evident. The structure in Fig. [TOl (left panel) 
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therefore does not imply significant correlation on the 
actual azimuth difference (/)a, no inherent collimation ef- 
fect in the initial state which might for example explain 
the SS ridge in p-p collisions at 7 TeV as suggested in 
Ref. [17]. In contrast, measured p-p hadron correlations 
depend strongly on 0a for hadron pt > 0.5 GeV/c and 
weakly for pt < 0.5 GeV/c (hard and soft components re- 
spectively), again contradicting the Glasma prediction. 



F. Glasma azimuth correlations and radial boosts 

A central issue for this study is the Glasma-model 
conjecture that the /^-elongated SS 2D peak observed in 
minimum-bias angular correlations ("soft ridge") repre- 
sents formation of a ID ridge on azimuth from boost 
of Glasma flux tubes emitting isotropically in their rest 
frames. Given the boost kinematics the peak on azimuth 
should be narrower for lower-p^ gluons (hadrons), and 
conversely for higher-p^ gluons (hadrons). That trend is 
opposite to what is actually observed for dihadron num- 
ber correlations with applied pt cuts [1, [29| , and for pt 
angular correlations j30|, [3l|. In contrast, the observed 
data trends are well explained by jet formation, where 
higher-pt particles contribute a narrower SS peak struc- 
ture. Boosted Glasma flux tubes are thus contradicted by 
Pt systematics of the SS 2D peak width. Furthermore, it 
was shown in Ref. [38| that a very large and problematic 
radial boost velocity {pt = 0.96) is required to achieve 
the narrow azimuth width observed in the data. 

The detailed SS peak shape on azimuth also provides 
important evidence. The Glasma model requires averag- 
ing SS peak widths over a broad radial-boost distribution 
and a pt spectrum. Such averaging invariably leads to 
long tails relative to the fundamental peak distribution 
(e.g. Sec. IV- A of Ref. [39j). No such tails are observed in 
the data. A narrow ID peak on azimuth with near-ideal 
Gaussian shape is observed for all Au- Au centralities [6] . 



G. Glasma gluon pseudorapidity correlations 

In Fig. 4 (left panel) of Ref. [l7| a Glasma prediction 
for gluon correlations on r]A is compared with measured 
Au-Au hadron correlations [40|. The Glasma prediction 
is essentially uniform on t^a for t^a ^ 4]. A "short- 
range" correlation peak introduced from PYTHIA {p-p 
collisions) to accommodate the heavy ion data is extra- 
neous to the Glasma prediction. No explanation is given 
for how a Glasma-based prediction of gluon correlations 
is related quantitatively to triggered dihadron correla- 
tions. Predicted variations at larger t^a are outside the 
T] acceptance of most detectors. 

Measured correlation structure on t^a provides essen- 
tial model tests. For all minimum-bias (p^-integral) an- 
gular correlations from all Au-Au centralities measured 
within the STAR TPC acceptance the SS 2D peak is 
consistent with a 2D Gaussian with large curvature on 



T^A (e.g. [Sl, l9|). In the static Glasma flux-tube scenario 
there is no variation on rj. The measured SS curvature on 
77a niust then (in the Glasma model) result from strong 
T] dependence of conjectured radial flow which is said to 
drive appearance of the SS 2D peak. However, the nonjet 
quadrupole nominally associated with elliptic flow shows 
no such T] dependence within the STAR TPC accep- 
tance [2l|. If a hydro interpretation is imposed on both 
the SS 2D peak (radial flow) and the nonjet quadrupole 
(elliptic flow) a major discrepancy emerges between two 
hydro manifestations. The Glasma model and jet-related 
T^A structure are further discussed in Ref. jl6| . 



H. Glasma theory comparison with hadron data 

In Eq. 5.3 of Ref. JlSj angular correlation data in the 
form Ap/y^p^ j6|, Q are compared with Glasma flux- 
tube predictions in the form 



^/Pref 



(0A = 0) 



dN Ap 

dy Pref 
1^2 



7b 



1 

IB 

13.5 as V^^ IB 



(6) 



where the LHS is evaluated as l/y27rcr^^ = 0.62 with 

=0.64. But that is the amplitude for a unit-normal 
Gaussian, not what was actually measured in Refs. [6|,[9[. 
The boost factor including is applied without justi- 
fication to correlation amplitude Ap/pref which has no 
structure on ^a- A SS 2D peak with large curvature on 
7^ A is implicitly compared to a ID ridge uniform on t^a- 
And the Glasma prediction scales as Npart whereas the 
measured peak amplitude increases faster than N^^i^. 

Since Ap/pref = /^2/A^FT, Eq. (|6]) as written implies 
that the gluon density per flux tube is {l/Npr) dNg/dy = 
1/13.5^5 ~ 1/7, since it is assumed that as ~ 0.5 {Q = 
0.8 GeV [41]). However, if factor l/27r which belongs in 
the first line (to match the LHS definition in Sec. Ill Cp 
is restored we obtain {l/Npr) dNg/dy ^ 1 compared to 
hadron duch/dr] = 2.5 for NSD p-p collisions (Lund string 
fragmentation). If we now insert those values into the 
first line of Eq. (|6]) we obtain 



(SSpeak) = 



or K2 



1 dNg 
Nft dy 27r 
0.7 X 27r 



IB 



1 

7b 



(7) 



7b 



1 

IB 



The estimate is larger by factor 2tt than that from 
Ref. [18|. Aside from the comparison method and the as- 
sumed radial flow boost the Glasma prediction is a factor 
27r too small relative to measured angular correlations, 
whatever the Au-Au centrality (unspecified in Ref. ^18j). 
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VI. DISCUSSION 

According to Ref. [17} correlations best reveal the A-A 
initial state (IS) dynamics (compared to integral yields 
and spectra). The Glasma picture of the IS in heavy ion 
collisions should be the natural framework to understand 
correlation mechanisms. In particular, Glasma flux tubes 
should explain the "ridge" encountered in more-central 
Au-Au collisions. To test that conjecture we contrast 
Glasma theory predictions of gluon correlations from 
Refs. [13, [3 with alternative models and with hadron 
fluctuation and correlation data. 



A. Glasma and "bulk" hadron production 

The CGC model describes "bulk" hadron production 
in terms of a longitudinal color-field system (Glasma flux 
tubes) approximating a dense gluonic system near mid- 
rapidity in A-A collisions. There is no transverse parton 
dynamics in the model, no large-angle parton scattering 
and fragmentation to jets. The Glasma model is a one- 
component (soft) model. Mid-rapidity gluons are pro- 
duced by emission from independent Bose-Einstein radia- 
tors (flux tubes) longitudinally boost invariant over some 
interval and with transverse correlation length ^ l/Qs- 

The Glasma model is formally similar to the Lund 
string model . The probability distribution of large-x 
color-charge field sources is analogous to a parton distri- 
bution function (PDF) [l7|. However, the Glasma model 
does not provide an absolute prediction for hadron pro- 
duction, only relative trends on A-A centrality and en- 
ergy. The relative centrality trend (Sec. IVIEp is contra- 
dicted by spectrum data [16]. 

In contrast, the Lund string model does provide quan- 
titative predictions of several aspects of nonperturbative 
soft parti cle p roduction in p-p and more-peripheral A-A 
collisions [Isl, |lj • The soft component observed in NSD 
p-p collisions [Ij (where "bulk" particle production is un- 
likely) appears to play a role even in central Au-Au col- 
lisions [3| , describing spectrum data quantitatively when 
supplemented by a hard component (jet fragments) de- 
scribed by pQCD [8]. In more-central Au-Au collisions 
the hard component contributes about one third of the 
total hadron production [5]. Correlations in all colli- 
sion systems require a two-component model including 
both soft production (Lund strings ^ flux tubes) and 
(semi) hard parton scattering and fragmentation ^6j. 



B. Gluon correlations inferred from hadron data 

Measured hadron correlations near mid-rapidity 
(Sec. IIV Ap exhibit apparent local charge and momen- 
tum conservation consistent with single-gluon parents 
common to two or more daughter hadrons, favoring a 
(1 — > 2) process. Hadron correlations seem to "point 



back" to single-parton (gluon) precursors. Hadron cor- 
relations from (2^1) coalescence would require that 
several parent partons conspire to produce observed two- 
hadron local net-charge and momentum correlations. 

Hadron spectrum structure can be combined with an- 
gular and pt X pt correlations consistent with the two- 
component model of hadron production to reconstruct 
the parent parton population near mid-rapidity. Hadron 
correlation data suggest that partons from dissociated 
projectile nucleons comprise a minimally-correlated low- 
Pt gluon population fragmenting longitudinally to charge- 
neutral hadron pairs which locally also conserve trans- 
verse momentum. That population produces the hadron 
correlation soft component which indicates no significant 
correlation among parent gluons. 

A scattered-gluon spectrum near mid-rapidity pre- 
dicted quantitatively by pQCD [H can be combined with 
measured fragmentation functions [JT] to predict hadron 
fragment distributions consistent with measured hadron 
hard components Q , thus confirming a parton spectrum 
with lower bound near 3 GeV. Hadron correlation struc- 
ture on Pt X Pt or yt x yt is consistent with that un- 
derlying parton spectrum. Analysis of hadron angular 
correlations indicates that scattered partons are corre- 
lated as momentum-conserving back-to-back recoil pairs, 
including an acoplanarity distribution depending on the 
initial-state parton kt spectrum, but are otherwise uncor- 
related. Even in central Au-Au collisions the pt spectrum 
of SS 2D peak hadrons and correspondence with the AS 
ID ridge compel interpreting the SS 2D peak in terms of 
large- angle scattering of energetic (> 3 GeV) gluons. 



C. Glasma gluon correlations on pt or yt 

Hadron correlations on pt x pt or yt x yt impose key 
constraints on any theoretical attempt to describe the 
SS 2D peak structure, yt x yt correlations from p-p col- 
lisions have a detailed quantitative correspondence to 
pQCD calculations and to spectrum hard components 
(Sec. IIV Ap . Decomposition of yt x yt correlation struc- 
ture according to combinations of SS and AS azimuth 
subregions and LS and US charge combinations reveals 
quantitative correspondence with expected features of 
parton scattering and fragmentation p, [sl, 3, 0, H, Ell . 

Per-pair measure K2 from Ref. [17] plotted in Fig. [3 
already reveals major disagreement between Glasma glu- 
ons and hadron per-pair p-p correlations from Ref. |3^ . 
Whereas the latter greatly increase in amplitude with 
larger yt the Glasma prediction in Fig. [3 generally de- 
creases to a constant asymptotic value near unity. If both 
results are converted to per-particle measures as in Fig. [8] 
the large discrepancy is again apparent, but the mea- 
sured hard-component structure appearing in the right 
panel [1, 01 corresponds quantitatively with pt spectrum 
structure jl|, [sj and pQCD calculations [8[. 

Because it lacks a parton large-angle scattering mecha- 
nism the Glasma model cannot describe measured yt x yt 
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correlations in p-p or more-central Au-Au collisions. The 
Au-Au case is quantitatively different from p-p colli- 
sions but still qualitatively incompatible with the Glasma 
model. Ironically, the Glasma model produces a very 
hard gluon pt spectrum as in Eq. ([3]), but nothing that 
corresponds to the observed hadron spectrum hard com- 
ponent which scales with A/'^^^, not Npart- 

D. Glasma gluon angular correlations 

Figure [10] (left panel) gives the impression that static- 
model gluon correlations are sharply peaked near paral- 
lel and antiparallel momentum configurations, suggesting 
that a SS ridge ( "collimation" ) is already inherent in the 
Glasma model without invoking radial fiow. In Sec. IV El 
the apparent peaked structure is identified as the Jaco- 
bian for the transformation 0a cos(0a)- Glasma cor- 
relations on 7^ A are also structureless within |7^a| < 4. 
Thus, the static Glasma model predicts no significant 
gluon angular correlations for any conditions, in contrast 
to measured p-p and Au-Au hadron angular correlations. 

According to the theory any hadron correlation struc- 
ture on azimuth must result from a conjectured radial 
boost. However, formation of a SS ridge via radial fiow is 
inconsistent with hadron data. Coupling boost-invariant 
fiux tubes with boost-invariant radial fiow would produce 
correlations on 0a alone. The measured large SS curva- 
tures on T^A would require strong variation of radial fiow 
on contradicting the measured uniform nonjet ^zimuth 
quadrupole within the STAR TPC acceptance j21]. 

The radial boost mechanism should produce a nar- 
rower structure on azimuth for lower-p^ particles. The 
measured hadron SS peak is narrower on azimuth for 
higher-pt particles [4], consistent with expectations for jet 
formation. There is no mechanism in the Glasma model 
for formation of the AS ID peak on azimuth which tracks 
very closely with SS 2D peak properties vs A- A central- 
ity. In contrast, the AS ID peak and its systematics are 
consistent with expected back-to-back jet correlations. 



Thus, centrality dependence of the SS 2D peak amplitude 
is incompatible with Glasma expectations. 

In the Glasma model the observed SS peak structure 
narrow on azimuth must result from conjectured radial 
fiow (Sec. IV Ep . But the magnitude of the radial boost 
and its effect on fiux tube emission must conspire to pro- 
duce the same SS 2D peak azimuth width ^ 0.65 
over a broad range of Au-Au centralities where the re- 
ported radial fiow magnitude Pt is changing from 0.25 to 
0.6 [36]. And the azimuth width must remain constant 
through and beyond the sharp transition where the am- 
plitude and T] width of the SS 2D peak (which must also 
depend on radial fiow) change rapidly [6]. 

The confiicting constraints on conjectured radial 
fiow systematics are inconsistent with observed nonjet 
quadrupole V2{2D} systematics which show no corre- 
spondence to the sharp transition in SS 2D peak char- 
acteristics [21]. The nonjet quadrupole is uniform on 
?7a within the STAR TPC acceptance. Any radial fiow 
would have to be very nonuniform within the same ac- 
ceptance to produce the SS 2D peak curvature on t^a- If 
the two phenomena are hydro manifestations they seem 
to be incompatible. A detailed search for radial fiow with 
two-component spectrum analysis was unsuccessful [3,]. 

F. Glasma energy dependence 

Per-particle fiuctuation measure {aj^ — N)/N = N/k 
represents an integral of angular correlations [20]. Thus, 
NBD parameter k in the form N/k represents a cor- 
relation integral of differential per-particle number an- 
gular correlations, including the SS 2D peak nominally 
described by Glasma fiux tubes. The predicted energy 
trend for parameter k in the Glasma model is /c ~ Nft oc 
^/s^ increasing monotonically with ^/s, implying that 
1/k should decrease slightly with energy. The observed 
trend in Fig. [6] is 1/k ex [ln(v/s/13.5 GeV)]^•^^ and 
the SS 2D peak amplitude itself is observed to increase 
with energy approximately as In (y^/ 13. 5 GeV) @, I3Q| . 
strongly contradicting Glasma expectations. 



E. Glasma centrality dependence 

The centrality dependence of measured spectra and 
correlations from Au-Au collisions has been extensively 
studied [1 H, i, §, [33|. The systematics of three SS 2D 
peak parameters challenge the Glasma fiux-tube model. 
The SS 2D peak (hard component) amplitude scales 
as Nijin in more-peripheral A-A collisions (as predicted 
for jets) and increases more rapidly than N^in in more- 
central Au-Au collisions above a sharp transition on cen- 
trality [6]. The Glasma fiux-tube number Nft scales 



with A-A centrality as QiS± 



N, 



part' 



Any Glasma 



fiux-tube contribution to the per-hadron SS 2D peak 
should actually decrease with increasing A-A centrality 
as Npart/'^'^ch ^ 1/[1 -\- x{u — 1)]. The measured SS peak 
amplitude increases at least as fast as h'/[l-\-x{u — l)] Q. 



G. Glasma gluon fluctuations and correlations 

The statistics of generic two-tiered particle production 
can be described as follows: If N particles result from 
event- wise production of K sources, each emitting n par- 
ticles (all are fiuctuating random variables) then [IF 



N 



N 



cr:: 



N 



+ A- 



K 



K 



-n,(8) 



illustratin g th e additivity of per-particle (co) variance 
measures [20|. It is assumed that K-n covariance is 
zero. That expression could describe gluon emission from 
Glasma fiux tubes, hadron emission from Lund strings 
or fragmentation of large-angle-scattered partons. The 
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terms on the RHS are interpreted in terms of jet pro- 
duction in Ref. [16] (Sec. V-A). The first term represents 
mtra source correlations, the second m^er source corre- 
lations, the third represents non-Poisson source number 
fluctuations and the last term represents the hierarchy 
process apart from any source or particle correlations. 

In the Glasma model the following are equivalent: 
K ^ k ^ Nft ~ Npart the number of flux tubes, 
n ^ Ug^ the number of gluons per independent color 
source, and N ^ Ng ^ rigNpr the total radiated gluon 
number. Independent color sources flux tubes) are 
assumed to be Poisson distributed with no intersource 
correlations, in which case the second and third terms on 
the RHS of Eq. ([8]) are zero. We then have from Eq. (j2j) 



f^2 



Nft— Nft ' 



cr:: 



Pref 
fin 



9 



(9) 



where the first term in the second line is the per-pair 
measure of all intrasource (per flux tube) correlations, 
explaining why 1^2 tends to unity at larger p^: intrasource 
correlations go to zero in the perturbative limit. In the 
limit where exactly one gluon proceeds from one flux tube 
(no statistical hierarchy) = 0, = 1 and ^ 0. 

Equation (|9]) contrasts dramatically with hadron data. 
Nft ^ Nparu ^2 = [Npart/ pQih)][pQ{h) /k] and both fac- 
tors have been measured. The first factor (centrality de- 
pendence) decreases ex 1/[1 + 0.1(i^ — 1)] [5]. The second 
factor (energy dependence) plotted in Fig. [6] (right panel) 
increases strongly with ^/s and is identified with several 
contributions from jet correlations in [I61] . 



H. Glasma flux tubes vs minijets in p-p and Au-Au 

The mechanism for ?7-elongation of the SS 2D peak in 
more-central Au-Au collisions [6, ''oj is a major problem 
for QCD theory at RHIC. Does the elongated peak arise 
from nonperturbative modification of parton scattering 
and fragmentation in a large A-A system, or does a novel 
process based on Glasma flux tubes and conjectured ra- 
dial flow manifest as an elongated ridge? Survival of 
copious parton scattering and fragmentation (minijets) 
in more-central Au-Au collisions [5|, l8|, l39|, |42j contradicts 
claims for formation of a flowing partonic medium with 
small viscosity [isl, Thus, an alternative explanation 
for the SS 2D peak by a mechanism other than parton 
scattering and fragmentation is sought. 

The SS 2D peak has been characterized as a "soft 
ridge" [32^ , and r] elongation is described as "long-range" 
correlations, confusing polar angle measure 77 and longi- 
tudinal momentum Pz represented by longitudinal rapid- 
ity yz- A causal argument is then invoked that only a 
process occuring at early times can produce long-range 
rapidity correlations. Glasma flux tubes are said to be in- 
herently long-range (boost invariant) and established at 



early times. They are thus characterized as a "natural" 
explanation for the 77-elongated SS 2D peak [l7j. 

However, while the Glasma model seems to describe 
qualitatively a few features of the SS 2D peak in more- 
central Au-Au collisions, there are substantial discrepan- 
cies between the flux-tube model and measured hadron 
spectrum and correlation systematics. Major issues in- 
clude: (non) existence of conjectured radial flow and in- 
consistency with observed azimuth quadrupole systemat- 
ics, disagreement with the energy and centrality depen- 
dence of hadron spectra and correlations, disagreement 
with most features of measured hadron pt and angular 
correlations, especially r] dependence of the SS 2D peak, 
absence of the pt x pt hard-component structure identi- 
fied with the 77-elongated SS 2D peak and absence of an 
AS ID ridge in the Glasma model. 

In contrast, the two-component soft+hard (=minijet) 
model of hadron spectra and correlations has been quan- 
titatively applied in combination with pQCD predic- 
tions to data from p-p and Au-Au collisions with good 
success [l|, @4El5 HI- The model continues to describe 
the SS 2D peak volume (sum of fragment pairs) for all 
Au-Au centralities. Large-angle parton scattering does 
occur at early times in the collision, with possible large- 
T] consequence in more-central A-A collisions. The only 
manifestation not currently accommodated by the two- 
component (fragmentation) model is rj elongation of the 
SS peak in more-central A-A collisions. However, devi- 
ations from pQCD expectations (e.g. angle asymmetries 
about the jet axis) are not inconsistent with conventional 
jet production in elementary collisions, for instance p-p 
collisions [7] and three-jet events in e+-e~ collisions [45|. 



VII. SUMMARY 

Observed strong elongation on 77 of the same-side 2D 
peak in minimum-bias angular correlations from Au-Au 
collisions has been attributed to Glasma flux tubes cou- 
pled with radial flow to form a narrow structure or ridge 
on azimuth. In the present study we have tested that 
conjecture by comparing Glasma predictions for particle 
production, spectra and correlations with conventional 
fragmentation models and with measurements. We find 
a number of contradictions between the Glasma model 
and spectrum and correlation data. 

The A-A centrality dependence of the Glasma model is 
defined by flux-tube number Nft^ approximated by the 
number of nucleon participants Npart- The Glasma one- 
component model of gluon (hadron) production therefore 
has no correspondence to observed hard-component fea- 
tures in spectra and correlations which scale with central- 
ity as the number of N-N binary collisions A/'^^^. There is 
no relation to p-p and peripheral A-A spectrum and cor- 
relation systematics described quantitatively by pQCD. 

The same-side 2D peak in angular correlations is asso- 
ciated with a prominent peak on pt x pt or yt x yt cor- 
relations with quantitative correspondence to the hard 
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component in measured pt spectra, all scaling as N^in. 
The hard-component peak remains visible near pt = 1 
GeV/c from p-p to central Au-Au collisions. The Glasma 
model predicts structure similar to observed soft compo- 
nent correlations which are unrelated to the same-side 
2D peak. No hard-component structure on pt x pt or 
yt X yt correlations is observed from the model, and the 
Glasma single-gluon spectrum scaling with Npart has no 
correspondence in hadron pt spectra. 

The static Glasma model exhibits no significant angu- 
lar correlations, relying on coupling to conjectured ra- 
dial flow to develop a peak structure on azimuth. How- 
ever, such a correlation mechanism would produce peak 
azimuth- width dependence on particle pt (e.g. trigger- 
associated cuts) opposite to that observed (which is con- 
sistent with pQCD jet structure). Radial flow would also 
have to exhibit strong variation on z to produce the large 
observed same-side peak curvature on t^a, which would 
be inconsistent with measured nonjet quadrupole system- 



atics. And the Glasma model has no mechanism to ac- 
count for the away-side ID peak on azimuth naturally ex- 
plained by back-to-back jet correlations scaling as N^n. 

The Glasma model does provide a qualitative conjec- 
ture as to a possible mechanism for same-side peak r] 
elongation. But Glasma predictions are contradicted by 
measured hadron spectra and correlations. In contrast, 
a two-component model of hadron production, includ- 
ing minimum-bias parton fragmentation within a pQCD 
framework, quantitatively describes most hadron data. 
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